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ABSTRACT 

The near-infrared to optical-ultraviolet (0.1 - 10 fj,m) spectral energy distribution 
(SED) shapes of 407 X-ray-selected radio-quiet type 1 AGN in the wide-field "Cosmic 
Evolution Survey" (COSMOS) have been studied for signs of evolution. For a sub- 
sample of 200 radio-quiet quasars with black hole mass estimation and host galaxy 
correction, we study the mean SEDs as a function of a broad range of redshift, bolo- 
mctric luminosity, black hole mass and Eddington ratio, and compare them with the 
Elvis et al. (1994, E94) type 1 AGN mean SED. The mean SEDs in each bin are very 
similar to each other, showing no evidence of dependence on any of the analyzed pa- 
rameters. We also checked the SED dispersion as a function of these four parameters, 
and found no significant dependance. The dispersion of the XMM-COSMOS SEDs is 
generally larger than E94 SED dispersion in the ultraviolet, which might be largely 
due to the broader "window function" for COSMOS quasars, and the X-ray based 
selection technique. 

Key words: galaxies: evolution; quasars: general; surveys 



1 INTRODUCTION 

The nature of the continuum emission of Active Galactic 
Nuclei (AGN) remains unsettled after several decades of 
study. Yet, understanding the continuum emission of AGN, 
* E-mail:hcnghao@post. harvard.edu from X-rays to radio, is essential to understand the physics 
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of accretion onto super massive black holes (SMBHs). The 
continuum in each spectral region can be ascribed to dis- 
tinct energy generation mechanisms: jets in the radio (see 
e.g. the review by Harris & Krawczynski 2006), dust in the 
infrared (IR, McAlary & Rieke 1988; Sanders et al. 1989), 
accretion disks in the optical-ultraviolet (UV) and soft X- 
rays (Shakura & Sunyaev 1973; Rees 1984; Czerny & Elvis 
1987), and Compton up-scattering by hot coronae in the 
hard X-rays (e.g., Zamorani et al. 1981; Laor et al. 1990; 
Haardt & Maraschi 1991; Williams et al. 1992; Zdziarski et 
al. 2000; Kawaguchi et al. 2001; Mateos et al. 2005; Mainieri 
et al. 2007). 

Continuum changes with redshift, luminosity or Ed- 
dington ratio might be expected. Most SMBH growth 
occurs during the active 'AGN' phases (the 'Soltan ar- 
gument', Soltan 1982), implying that most bulges went 
through an AGN phase (e.g. Magorrian et al. 1998). Rapid 
growth of central SMBHs happens in high-redshift and high- 
luminosity quasars emitting near the Eddington limit (Barth 
et al. 2003; Vestergaard 2004; Jiang et al. 2006; Kurk et al. 
2007; but see also Steinhardt & Elvis 2010). The space den- 
sity of X-ray-selected, highly luminous AGN peaks at around 
z = 2.5, and declines at z > 3 (Silverman et al. 2005, Brusa 
et al. 2009, Civano et al. 2011). Low- luminosity AGN are 
more prevalent at z < 1 than higher luminosity ones (Cowie 
et al. 2003; Fiorc et al. 2003; Ueda et al. 2003; Silverman et 
al. 2005) . As the central SMBH is the driver of the emission, 
one might expect the quasar SED to evolve as the black hole 
grows due to accretion. 

Many parameters - the black hole mass, the AGN lu- 
minosity relative to the host galaxy, the accretion rate, the 
physical properties of the accretion disk and the properties 
of the absorbing dust - might affect the shape of the AGN 
SED (Wilkes 2003). For instance, the optical to X-ray spec- 
tral index [ctox =0.3841og(F2fc e v /F °)], correlates with 
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luminosity but not with redshift (e.g. Vignali et al. 2003a; 
Steffen et al. 2006; Just et al. 2007; Young et al. 2010; Lusso 
et al. 2010). It is possible that similar dependency of lumi- 
nosity in the SED shape exists at other wavelengths. 

Observations indicate a tight link between SMBH 
growth and galaxy evolution (e.g. Magorrian et al. 1998; 
Marconi & Hunt 2003; Tremaine et al. 2002; Menci et al. 
2008). Locally, SMBHs appear to reside at the center of 
most galaxies and the SMBH masses are tightly correlated 
with their masses (e.g. Kormendy & Richstone 1995; Mar- 
coni & Hunt 2003) and velocity dispersions (i.e. Mbh — a 
relations; Ferrarese & Merrit 2000; Gebhardt et al. 2000; 
Tremaine et al. 2002). Some evidence for evolution of this 
relationship has been reported using several methods (e.g. 
Peng et al. 2006; Shields et al. 2006; Ho 2007; Merloni et al. 
2010). This evolution would imply that the feedback of the 
SMBH to the host galaxy evolves. Similar evolution in the 
innermost regions (within the 'torus') is possible too. All of 
these processes could lead to different SED shapes. 

However, no evolution of the AGN SED has yet been 
demonstrated. There is no convincing evidence for any 
change of SED with redshift (Silverman et al. 2002; Mathur 
et al. 2002; Brandt et al. 2002; Vignali et al. 2003b). High 
redshift quasars (up to redshift 7) show optical spectra simi- 
lar to low redshift quasars from the SDSS (Jiang et al. 2007, 
Mortlock et al. 2011). There is evidence, though, that the 
SEDs of extremely low-luminosity (Lboi < 10 42 erg/ s) nu- 



clei are remarkably different from those of luminous (Lboi > 
10 44 erg I a) AGN (Ho 1999, Ho 2008). 

So far, the systematic study of the dependency of the 
SED shape with physical parameters has been limited by 
difficulty in obtaining a large sample size with good multi- 
wavelength coverage. The Cosmological Evolution Survey 
(COSMOS, Scoville et al. 2007) has the appropriate combi- 
nation of depth, area and extensive multi-wavelength data 
that allows for a sensitive survey of AGN to address this 
question. 

The COSMOS field has been imaged with XMM- 
Newton for a total of ~ 1.5 Ms (Hasinger et al. 2007; Cap- 
pelluti et al. 2007, 2009). Optical identifications were made 
by Brusa et al. (2010) for the entire XMM-COSMOS sam- 
ple. Photometric properties and redshifts were produced for 
each point source. This extensive data set allows us to make 
a systematic study of the evolution of the SED shape, which 
is the main purpose of this paper. 

From this complete sample, we extracted a sample of 
413 type 1 AGN (broad emission line FWHM>2000 km s _1 ). 
The type 1 AGN SED sample catalog is described in detail 
in Elvis et al. (2012, hereafter Paper I). It includes quasars 
with redshifts 0.1 ^ z ^ 4.3 and magnitudes 16.9 ^ iab ^ 
24.8, with 98% of the sources being radio-quiet (Hao et al. 
2012a). This sample is twenty times larger than the Elvis 
et al. (1994, E94 hereafter) radio-quiet type 1 AGN SED 
sample, and has full wavelength coverage from radio to X- 
rays (for a total of 43 photometric bands, Paper I) and high 
confidence level spectroscopic redshifts (Trump et al. 2009a; 
Schneider et al. 2007; Lilly et al. 2007, 2009). The mean 
SED of the XMM-COSMOS type 1 AGN was calculated and 
compared to previous studies: E94, Richards et al. (2006), 
Hopkins et al. (2007) and Shang et al. (2011) in paper I. We 
found that in the near-IR to optical-UV bands, the previous 
studies have very similar shapes, while the XMM-COSMOS 
mean host-corrected quasar SED has a less prominent 'big- 
blue bump', due to possible excess host contribution, not 
corrected because of the dispersion in the black hole mass 
and host luminosity scaling relationship itself. In this paper, 
we compare the Paper I sample with E94 as a representative. 

Paper I presented the selection and properties of 
the XMM-COSMOS type 1 AGN sample of 413 quasars 
(XC413 hereinafter). We used seven radio- loud criteria 
(Rl = log(/ 5 GW/s) > 1. Wilkes & Elvis 1987; q 2 4 = 
log(/24 M m//i.4GHz) < 0, Appleton et al. 2004; R uv = 
log(/ 5G W/ 2500 ^) > 1, Stocke et al. 1992; P 5G h z = 
\og[P 5G Hz(W/Hz/Sr)] > 24, Goldschmidt et al. 1999; and 
Rx = \og(vL v (5GHz)/L x ) > -3, Terashima & Wilson 
2003) to define a radio-loud quasar. We find that the radio- 
loud fraction is 1.5%-4.5% using any criterion, except R uv , 
which is subject to the reddening and host contamination 
issues (Hao et al. 2012a). Using two criteria at the same 
time, the radio-loud fraction is < 8/413 = 2%. Only 6 
XC413 quasars satisfy all the seven criteria. We define these 
6 quasars in this catalog as radio-loud (Paper I; Hao et 
al. 2012a). We refer to the radio-quiet XC413 as XCRQ407 
hereinafter. 

Estimates of black hole mass (Mbh) for 206 of XC413 
have been made by Merloni et al. (2010) and Trump et al. 
(2009b). Both papers used single-epoch spectra and applied 
the scaling relations from reverberation mapping of Vester- 
gaard & Peterson (2006). This method requires a high S/N 
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Table 1. Source Properties in Parameter Space 



XID tP_ logMBH^ \ogL bo £_ logAj^ \ogL a i£_ L ir /L bol f_ L p t /L bol 9_ L x /L bo £_ logL bol:ht S_ 







[M Q ] 


[erg s 1 ] 




[erg s 1 ] 


% 


% 


% 


[erg s 


1 


0.373 


8.58 


45.45 


-1.22 


45.46 


19.7 


38.2 


5.8 




2 


1.024 


8.96 


46.02 


-1.04 


46.09 


38.2 


20.9 


32.4 


46.00 


3 


0.345 


8.66 


45.37 


-1.39 


45.39 


35.6 


25.4 


18.2 


45.28 


■1 


0.132 


7.31 


44.50 


-0.91 


44.55 


45.7 


39.7 


3.4 


44.43 


5 


1.157 




46.16 




46.21 


38.5 


28.5 


14.5 




6 


0.360 


8.64 


45.10 


-1.64 


45.20 


39.2 


43.7 


6.0 


44.91 


7 


0.519 


8.38 


45.42 


-1.06 


45.44 


39.5 


28.9 


18.8 


45.37 


8 


0.699 


7.96 


45.96 


-0.11 


45.97 


34.1 


33.1 


7.5 


45.96 


9 


1.459 


8.86 


46.19 


-0.77 


46.28 


48.9 


22.4 


14.3 


46.18 



A portion of the table is shown here for guidance. The complete table will be available online. 



a The spectroscopic redshifts are from Trump et al. (2009a), Schneider et al. (2007), Lilly et al. (2007, 2009) 
^ The black hole mass estimates are from Trump et al. (2009b) and Merloni et al. (2010) 
c Calculated by integrating the rest frame SED from 24 fim to 40 keV. 
d \ E = L bo i/L Edd , details in §[23] 

e Calculated by integrating the rest frame SED from 1.4 GHz to 40 keV. 
f Li r is calculated by integrating the rest frame SED from 24 (im to 1 fim. 
9 Lopt is calculated by integrating the rest frame SED from 1 fim to 912A. 

Lx is calculated by integrating the rest frame SED from 0.5 keV to 40 keV. 
1 Calculated by integrating the rest frame host-corrected SED from 24 fim to 40 keV for 203 quasars in the XMM-COSMOS sample. 



spectrum with the broad emission line not near the ends of 
the spectra. For the quasars with only zCOSMOS spectra, 
the black hole mass was estimated only for those with Mgll 
lines in the spectra (Merloni et al. 2010), using the calibra- 
tion of McLure & Jarvis (2002). 206 quasars out of XC413, 
that is 203 quasars out of XCRQ407 have black hole mass es- 
timates. Paper I estimated the host galaxy contribution for 
203 of these 206 quasars using the scaling relationship be- 
tween black hole mass and host luminosity (Marconi & Hunt 
2003) adding an evolutionary term (Bennert et al. 2010, 
2011), excluding the 3 which had over-subtraction problems 
(the estimated host galaxy luminosity is larger than the ob- 
served luminosity). We define this sub-sample (SS) as SS203. 
In SS203, 200 quasars are radio-quiet. We refer to this sub- 
sample as SSRQ200. 

In Paper I, rest frame SEDs of XC413 were constructed 
on a uniform frequency grid (A log!/ = 0.02) from radio 
to X-rays. As there is limited data in the far-IR (> 10 fim 
in rest frame) and longer wavelengths, and because the UV 
flux is strongly affected by both variability and strong broad 
emission lines, we defer the analysis of these regions to a 
later paper. In this paper, we focus on the SED shape in the 
optical to near-IR range, while all the plots are shown in 
the rest frame 0.1 — 10 fim range. Note that all the sources 
discussed in this paper are broad line AGN with typical 
luminosity (~ 10 44 — 10 47 erg/s), typical black hole mass (~ 
10 7 — 10 9 Mq) and typical accretion rates (Eddington ratio 
larger than 0.01, see in S! 12. 3[) . 

All the wavelengths discussed in this paper are in the 
rest frame. We adopt the WMAP 5-year cosmology (Ko- 
matsu et al. 2009), with H =71 km s" 1 Mpc" 1 , Q M = 0.26 
and « A = 0.74. 



Table 2. Parameter Range 



Parameter 


XCRQ407 


SSRQ200 




min 


med max 


min 


med max 


z 


0.10 


1.57 4.26 


0.13 


1.50 4.26 


iog(M BH /M Q ) 






7.18 


8.39 9.34 


log-Laii 


44.29 


45.71 47.40 


44.55 


45.77 47.40 


logifcoi 


44.26 


45.69 47.28 


44.50 


45.74 47.28 


Lir / L bo { 


12.8% 


35.1% 82.5% 


12.8% 


35.6% 82.5% 


Lopt 1 L bo i 


9.3% 


33.3% 61.1% 


9.3% 


33.3% 59.7% 


Lx/Lbol 


0.34% 


8.6% 45.0% 


0.48% 


8.4% 45.0% 


logA B 






-1.76 


-0.72 0.63 



2 PARAMETER SPACE 
2.1 Bolometric Luminosity 

The bolometric luminosity is the total energy radiated by 
the quasar at all wavelengths. This luminosity is, in princi- 
ple, simply calculated by directly integrating the rest-frame 
SED over the whole wavelength range. In practice, this is 
difficult and observationally expensive. For the sample dis- 
cussed in this paper, it is possible to approximate the bolo- 
metric luminosity with the SEDs available (Paper I). 

As described in Paper I, the SEDs are produced by lin- 
early interpolating between the data points in log vL v ver- 
sus log v space (i.e. connecting the individual points with 
power laws in linear space). The COSMOS photometry for 
the XC413 sample is >90% complete from u (CFHT) to 
MIPS 24/im, that is over the 1.8 dex wide 0.35/im - 24/im 
(~0.14mum - 10/im for the typical z = 1.5 of XC413) ob- 
served frame interpolation is unproblematic (Paper I). In 
the mid-infrared range, for quasars with 70/im or 160/im 
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detections, we joined the 24^im data to the longer wave- 
length points with a power-law in log vi v vs. log v space; 
for the others we extrapolated from the rest frame 24/im 
to 8/im slope and checked that the extrapolation generally 
works. As the far-IR to radio photometry data are sparse, 
for each source with a > 3<j VLA detection, we assumed a 
power law f v tx i/ -0 ' 5 (e.g., Ivezic 2004) in the rest frame 
1.4 GHz (21 cm) to 100 GHz (3 mm) range. In the 100 GHz 
(3 mm) to 160/xm part, the SED can be approximated by 
the red end of the grey body /„ tx v 3+0 /(e h " /kT - 1), when 
hu < kT, f v tx u 2+l3 (e.g. Lapi et al. 2011). The /3 value is 
generally chosen in the range 1-2 (Dunne & Eales 2001). We 
thus choose j3 — 1, i.e. we assume a power-law /„ tx v 3 in 
this wavelength range. We directly linearly interpolate from 
Lyman break (1216A) and 0.5 keV SED (Laor et al. 1997) 
in log uL v versus log v space. In the X-rays, we use the mea- 
sured photon spectral index (F, f v cx i/ 1_r ' ) Mainieri et al. 
2007) and the observed 2 keV luminosity to get the SED in 
the 0.5 keV-40 keV (rest-frame), which is the range we could 
have XMM data for XC413. 

We used the SEDs from Paper I to calculate two ap- 
proximations to the bolometric luminosities by integrating 
the rest frame SED over several different wavelength ranges: 
(1) from 24 fim to 40 keV, the range over which COS- 
MOS coverage is virtually complete, as Lboi and (2) from 
21 cm/1.4 GHz to 40 keV as L a ii- L a ii was integrated over 
all the wavelengths for which we could have data. 

As the great majority of the XC413 quasars are radio- 
quiet and the photometric coverage from 24 \im to 1.4 GHz 
is currently sparse (Paper I), we use Lboi to avoid unwar- 
ranted extrapolation. The contribution of the radio emis- 
sion to the bolometric luminosity is less than 3% even for 
the radio-loud sources. Most of the far-IR (at > 100/im) 
luminosity is probably due to star-forming activity not the 
AGN (Netzer et al. 2007, Mullaney et al. 2012; but see Ho 
2005). Hence Lboi is a good approximation of the true bolo- 
metric quasar luminosity defined over the entire wavelength 
range. 

We note that the Lboi calculated by directly integrating 
the SED in the rest frame from IR to X-rays overestimates 
the quasar emission for two reasons: (1) the contribution 
from the host galaxy is not excluded; (2) different wave- 
length ranges of the SED have different emission origins. 
From the optical short-ward the primary emission is from 
the innermost region of the quasar (SMBH and accretion 
disk). The near and mid-IR continuum is the result of repro- 
cessing shorter wavelength radiation of the quasar by dust 
(Sanders et al. 1989; Suganuma et al. 2006). The reprocessed 
dust emission thereby includes reprocessed primary radia- 
tion emitted in directions different from our line-of-sight. In 
a non-spherical geometry, as is likely present in quasars, this 
makes us count more radiation that we would in a 4-7T aver- 
aged calculation. In this sense, we are double-counting the 
primary emission (e.g. Lusso et al. 2010). 

In practice, we do not have good estimates of the host 
contribution. For SSRQ200, we calculated the bolometric 
luminosity by integrating the host corrected SED from 24/xm 
to 40 keV as Lboi,hc (last column in Table [T}. As discussed 
in Paper I, we applied the scaling relationship between black 
hole mass and host luminosity reported in Marconi & Hunt 
(2003), and added an evolutionary term (Bennert et al. 2010, 
2011), to estimate the contribution of host galaxy. 




13.5 14 14.5 15 15.5 

\ogv (Hz rest) 



Figure 1. The frequency dependent bolometric correction for 
the SSRQ200 sample (black solid). The bolometric corrections in 
different redshift bins are also shown in the same plot (redshift: 
0.1-1.2= magenta short dashed line, 50 quasars; 1.2-1.5= cyan 
long dashed line, 50 quasars; 1.5—1.8= blue dot dashed line, 51 
quasars; and 1.8-4.3= green short long dashed line, 49 quasars). 
The red solid line show the bolometric correction for E94 RQ 
mean SED. The red star point at B band (4400 A) show the bolo- 
metric correction calculate from the Hopkins et al. (2007), assum- 
ing logL(, ; = 45.5 erg/s, which is the median L(, ; of SSRQ200. 
The red square show the bolometric correction at 3000A and 
5100A respectively from Runnoe et al. (2012). In the bottom 
panel we show the dispersion of the bolometric correction for the 
SSRQ200 sample. 



As mentioned above Lboi might overestimate the quasar 
emission for two reasons, meanwhile Lboi might also under- 
estimate the total emission for three reasons: (1) excluding 
the far-IR to radio emission; (2) excluding the hard X-ray 
emission, where we do not have data; (3) ignoring the red- 
dening of the primary emission in the optical/UV. Therefore, 
in this paper, we necessarily use Lboi as an estimate of the 
emission from everything associated with quasar hereafter, 
including emission from the primary radiation as well as the 
radiation reprocessing of the primary emission. 

To investigate the contribution to Lboi by different 
wavelength ranges, we calculated Li r by integrating from 
24 fj,m to 1 )im, L op t from 1 fim to 912 A and Lx from 
0.5 keV to 40 keV. The fractions of the luminosity in dif- 
ferent wavelength ranges are reported in Table [1] (The full 
table is available on line). Table [2] shows the median values 
and the ranges for Lboi, the fractions in three wavelength 
ranges, and the median values and the ranges for Mbh and 
Eddington ratio (which will be discussed in § 12 . 3f) . In gen- 
eral, the IR component (1 — 24 /im) provides an equal or 
greater contribution compared to the optical/UV 'big blue 
bump' component (0.1 — 1 /im). 



SEDs of Type 1 AGN in XMM-COSMOS II 5 

Table 3. Bolometric Corrections for the SSRQ200 sample 1 



logA logy All z 0.1 < z < 1.2 1.2 < z < 1.5 1.5 < z < 1.8 1.8 < z < 4.3 

(fim) (Hz) BC a BC a BC a BC a BC ct 



1.0 


13.5 


9.78 


5.37 


9.89 


4.61 


9.81 


5.82 


8.60 


3.97 


10.86 


6.63 


0.9 


13.6 


9.20 


4.40 


10.14 


4.45 


9.65 


4.53 


8.18 


3.33 


8.85 


5.01 


0.8 


13.7 


9.25 


4.05 


10.53 


4.54 


9.64 


3.92 


8.39 


3.03 


8.44 


4.27 


0.7 


13.8 


9.39 


3.90 


10.62 


4.49 


9.74 


3.69 


8.67 


2.93 


8.53 


4.10 


0.6 


13.9 


9.55 


3.72 


10.23 


3.84 


10.03 


3.76 


9.10 


3.19 


8.83 


3.97 


0.5 


14.0 


9.32 


3.53 


9.27 


3.50 


9.38 


3.19 


9.42 


3.55 


9.21 


3.93 


0. 1 


14.1 


8.66 


3.09 


8.75 


3.38 


8.59 


2.82 


8.15 


2.22 


9.18 


3.76 


0.3 


14.2 


8.95 


3.84 


9.36 


5.60 


9.29 


3.29 


8.22 


2.38 


8.93 


3.37 


0.2 


14.3 


10.11 


5.08 


10.41 


8.09 


10.67 


4.16 


9.18 


2.89 


10.21 


3.57 


0.1 


14.4 


11.24 


5.52 


11.74 


8.47 


12.45 


5.38 


10.18 


3.21 


10.61 


3.08 


0.0 


14.5 


12.94 


11.18 


14.20 


16.50 


15.63 


13.95 


11.28 


3.64 


10.64 


3.24 


-0.1 


14.6 


12.42 


16.46 


15.82 


31.83 


12.37 


6.52 


10.75 


3.17 


10.72 


4.13 


-0.2 


14.7 


11.31 


14.54 


15.59 


28.17 


10.97 


4.24 


8.97 


2.68 


9.73 


3.57 


-0.3 


14.8 


9.91 


5.62 


11.98 


9.37 


10.50 


3.70 


8.56 


2.58 


8.59 


3.44 


-0.4 


14.9 


8.62 


3.00 


9.73 


2.91 


9.58 


3.29 


7.86 


2.54 


7.29 


2.50 


-0.5 


15.0 


8.13 


4.25 


9.77 


5.23 


8.96 


4.59 


7.41 


3.46 


6.35 


2.38 


-0.6 


15.1 


9.03 


9.44 


11.97 


16.46 


9.89 


6.80 


7.49 


3.92 


6.75 


3.52 


-0.7 


15.2 


9.64 


11.98 


13.63 


19.77 


10.65 


11.24 


7.74 


4.51 


6.51 


3.77 


-0.8 


15.3 


9.43 


13.56 


14.09 


21.72 


10.38 


13.49 


6.94 


5.60 


6.28 


4.61 


-0.9 


15.4 


10.17 


14.21 


14.62 


23.45 


10.84 


13.06 


8.11 


5.96 


7.08 


5.38 


-1.0 


15.5 


12.85 


15.19 


16.03 


23.17 


11.82 


12.93 


9.91 


6.59 


13.72 


13.07 



1 After host galaxy subtraction. 



2.2 Bolometric Correction 

We then calculated the bolometric correction at various fre- 
quencies. The bolometric correction [BC V = Lboi/vL w ) is 
the factor which transforms the luminosity in one band to 
the bolometric luminosity. As the host galaxy contribution 
is prominent for X-ray-selected quasars, we calculated the 
frequency-dependent BC in the rest frame 0.1-10/xm only 
for SSRQ200 after performing the host correction. The bolo- 
metric correction was calculated for a A log v — 0.02 grid. 
The bolometric luminosity used in the calculation is L 00 i,hc, 
which is the integration of host-corrected SEDs, listed as the 
last column in Table [T] 

The mean and dispersion of the BC U are listed in Table 
[3] and shown in Figure [T] where we also show the mean and 
dispersion of the BC V for quasars at different redshift bins. 
The mean BC V curves for different redshift bins are consis- 
tent with each other given the large dispersion at each wave- 
length (Figure [1} . The dispersion is the largest at around 
1 (im, where the host contribution is the highest and in the 
extreme UV part where quasar variability is likely to con- 
tribute significantly to the observed dispersion. 

For comparison, we plot the E94 mean SED bolomet- 
ric correction as red solid line in Figure [T] Hopkins et al. 
(2007) used a double power law to approximate the B band 
bolometric correction. 

T / r \ -° 37 / t \ -0 012 

Lbol a r,r I L bol \ , n nn / Lb ° l \ 

To compare with the bolometric correction we got, we use 
the median value of Lbol (10 45 ' 5 erg/s) in SSRQ200 to apply 
to the formula. The B band bolometric correction from Hop- 
kins et al. (2007) is thus 9.99 (red star in Figured]). More 
recently, Runnoe et al. (2012) study the SED of 63 bright 
quasars at low redshift and got linear bolometric correction 
of 4.2 ±0.1, 5.2 ±0.2 and 8.1 ±0.4 at 1450, 3000 and 5100A, 



respectively (red square in Figure [T]) . These results gener- 
ally agree with bolometric correction got from SSRQ200. 
Note that the X-ray-selected SSRQ200 generally do not have 
prominent 'big-blue-bump' as the optical-selected E94 or 
Runnoe et al. (2012) sample, thus the bolometric correc- 
tion in the optical part is generally larger than the Runnoe 
et al. (2012) number. 

2.3 Eddington Ratio 

For SS200, the accretion rate relative to the Eddington rate 
(the Eddington ratio, Xe) can be calculated given the bolo- 
metric luminosities derived in i; 12.11 i.e. 

^ _ Lboi _ Lboi _ Lboi 

E ~ TsTd ~ ^g°"> Mbh ~ 1-26 x 1O 38 (M SH /M ) 

The black hole mass [Mbh) and the corresponding log- 
arithm of the Eddington ratio (logAg) of the XC413 are 
listed in Table [T] (Full table available on-line) . The median 
and ranges of these parameters are reported in Table [5] 

2.4 Parameter Space Classification 

The XCRQ407 sample spans a wide range of z, Lboi, Mbh 
and Xe (Tabled Figure[5]). The redshift range is comparable 
to that of Spitzer-SDSS sample (Richards et al. 2006), and 
the luminosity, black hole mass and Eddington ratio ranges 
are comparable to that of PG quasars (Sikora et al. 2007). 
The size of the XCRQ407 sample is more than double that 
of previous samples, and is an order of magnitude larger 
than the E94 sample. This large sample size spanning a wide 
range of the parameter spaces is useful to understand where 
and how the properties of the SMBH affect the SED shape. 
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Table 4. Parameter Space Bins 



bin 


redshift 


N 




Lbol 


N 




M BH 


N 


logA s 


N 






XCRQ407 


SSRQ200 




XCRQ407 


SSRQ200 




SSRQ200 




SSRQ200 


1 


0.103-1.166 


107 


50 


44.26-45.44 


125 


50 


7.18-8.10 


50 


-1.828- -1.040 


50 


2 


1.166-1.483 


72 


50 


45.44-45.74 


95 


50 


8.10-8.39 


50 


-1.040- -0.728 


50 


3 


1.483-1.848 


84 


51 


45.74-46.10 


103 


50 


8.39-8.66 


50 


-0.728- -0.375 


50 


4 


1.848-4.256 


144 


49 


46.10-47.29 


84 


50 


8.66-9.34 


50 


-0.375- 0.631 


50 




Redshift Redshift Redshift 

Figure 2. The classification of z, Lioi, Mbh and logAg bins of the XMM-COSMOS type 1 AGN sample. The dashed lines show the 
borders of quartile bins: the z bins (0.1 ~ 1.2 ~ 1.5 ~ 1.8 ~ 4.3), L bol bins (44.3 ~ 45.4 ~ 45.7 ~ 46.1 ~ 47.3), Mbh bins (7.2 ~ 8.1 
~ 8.4 ~ 8.7 ~ 9.4) and logA E bins (-1.8 ~ -1.0 ~ -0.7 ~ -0.4 - 0.6). The black crosses show all the radio-quiet XMM-COSMOS type 1 
AGN. The red triangles show the 6 radio-loud quasars. The green squares show the sub-sample SS203. 



We must note that XCRQ407 is an X-ray-selected sam- 
ple, and also includes sources with large host galaxy con- 
tribution (Paper I). For this reason, we focus on the host- 
corrected SEDs of the SSRQ200. For comparison we also 
study the SEDs of the XCRQ407 before the host correction. 
The median values and ranges in z, Lbol, Mbh and logAe 
are similar for the two sub-samples (Table [5]). 

We can now check the SED shape dependence on each 
physical parameter by dividing the samples into quartiles of 
quasar z, L bo i, Mbh and Xe, as shown by the dashed red 
lines in Figure [2] The bin boundaries and number of sources 
in different bins are given in Table [4] With these divisions, 
different bins have a similar number of quasars to calculate 
the mean SED, and therefore the possible statistical differ- 
ences between different bins are minimized. 



3 MEAN QUASAR SED DEPENDENCY ON 
PHYSICAL PARAMETERS 

We can look for trends in the SEDs by checking the mean 
SED shape diversity in different bins. For each of the z, 
Lboi, Mbh and logAe bins the mean SED are calculated as 
in Paper I. Briefly we: (1). converted the flux densities at 
each frequency for each object to luminosity, using a ACDM 
Concordance Cosmology (Komatsu et al. 2009); (2). shifted 
them to the rest frame for each source; (3). corrected for 
the small Galactic extinction (E(B-V)~0.017); (4). limited 
the variability by restricting the photometry data in use to 
2004 - 2007; (5). corrected for broad emission line contri- 
butions, which can be significant in the intermediate width 



Subaru bands; (6). linearly interpolated the SED to a uni- 
form frequency grid (Alog^ = 0.02); (7). calculated the 
mean SED at each log v grid point. To avoid the SED shape 
being dominated by the few luminous quasars in each bin, 
we also calculated the mean of the SED after normalizing 
the SED at 1/im. 

The resulting mean SEDs in different z and Lboi bins 
for XCRQ407 are shown in Figure For ease of comparison, 
we also plot the E94 mean SED and a galaxy template from 
SWIRE (Polletta et al. 2007) normalized to the value of 
L* from the UKIDSS Ultra Deep Survey (Cirasuolo et al. 
2007, M* K = -23). The galaxy SED shown is an Elliptical 
galaxy with an age of 5 Gyr (E5 hereafter). Different galaxy 
templates have similar shapes at around 1/im, so here we 
just choose E5 as a representative case. 

The SEDs in our sample have much less pronounced 
1/tm inflection point than in E94. In Paper I we concluded 
that this shape is probably due to the host galaxy contribu- 
tion. It is clear from Figure [3] that for low z and low Lboi 
sources the galaxy component around 1/tm strongly affects 
the shape of the SEDs. Indeed at higher z and Lboi, the 1/im 
inflection becomes more obvious, as the quasar component 
becomes relatively stronger. However, even for high z and 
Lboi sources, the optical big blue bump is not as strong as in 
E94. This is probably because E94 considered a UV selected 
sample which picks out the bluest quasars. The mean SED 
shapes in the central z and Lboi bins are quite similar to the 
mean SEDs of all the quasars in XCRQ407. 

For the SSRQ200 sample, we calculate the mean host- 
corrected SED for each bin, as shown in Figure [4] Even with 



SEDs of Type 1 AGN in XMM-COSMOS II 7 




13.5 14 14 5 15 15 5 13.5 14 14 5 15 15 5 13.5 14 14 5 15 15.5 13.5 14 14 5 15 15.5 



Figure 3. The mean SED for XCRQ407 in bins of z, L^oi before and after normalization at 1/im compared to E94 mean radio-quiet 
SED (red solid line). The black lines show the mean SED for all the radio-quiet type 1 AGN in the XMM-COSMOS sample. A host 
galaxy template (an Elliptical (5 Gyr), E5, blue solid line) is normalized to L„ from UKIDSS Ultra Deep Survey (Cirasuolo et al. 2007). 




Figure 4. The host-corrected mean SED for SSRQ200 in bins of z, L ao i, Mbh and A^; compared to E94 mean radio-quiet SED (red solid 
line):top: before normalization; bottom: normalized at 1/im. The black solid lines show the host-corrected mean SED of all quasars in 
SSRQ200. A host galaxy template (an Elliptical (5 Gyr), blue solid line) is normalized to from UKIDSS Ultra Deep Survey (Cirasuolo 
et al. 2007). Note that, to show the difference of the mean SEDs in different bins clearly, wc expanded the y axis in all the plots of the 
bottom panel. 



the host-correction, the mean SEDs for the low z, Lboi and 
Mbh sources are still relatively fiat and different from E94. 
This may be due to insufficient host galaxy correction in 
some cases because of the scatter in the scaling relationship. 
The mean SEDs show that the top three bins of z, Lboi and 
Mbh have similar shapes, showing no sign of dependency on 
any of these parameters. The 1 pm-normalized mean SED 
shapes of three log Xe bins are also similar to each other 
except for the penultimate bin (—1.0 < log A_e < —0.7). The 
z bins have the tightest distribution compared to other pa- 
rameters (e.g., at 3pm and 3000A the difference of the mean 
SEDs in different bins is 0.1 dex and 0.15 dex respectively, 
smaller than corresponding values of the other parameters). 
The difference of the mean SEDs in bins of each parameter 
are less than a factor of 2 even at the infared/UV end. 

To check for partial dependencies of the SED shape on 



the physical parameters, we checked the SED shape differ- 
ence with one parameter when fixing another. In order to 
have a large enough number of quasars in the mean SED 
calculation in each sub-bin, we consider the central two bins 
(Figure [5]). The number of quasars in each sub-bin is listed 
in Table [5] As the number of quasars in each sub-bin is 
relatively small, these mean SEDs are more affected by par- 
ticular SED shapes. In all these plots, the bottom / first bins 
are still affected by the host contribution, and the top / last 
bins are sometimes affected strongly by several specific SEDs 
(e.g. hot-dust-poor quasars, Hao et al. 2010, 2011). However, 
in general, the mean SEDs shapes in different sub-bins are 
very similar to each other and to the mean SED in that bin. 
For instance, when fixing z, the top / last two Lboi sub- 
bins have a similar mean SED shape, while the mean SED 
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Figure 5. The mean host-subtracted SED normalized at lfim compared to E94 mean radio-quiet SED (red solid line). Different rows 
are for quasars in specific bins: (f) z bin 2 (left half 1.2 < z < 1.5) and bin 3 (right half, 1.5 < 2 < 1.8); (2) logLi, ; bin 2 (left half 45.4 < 
logL,, oi < 45.7) and bin 3 (right half, 45.7 < logL boi < 46.1); (3) \og{M BH /M Q ) bin 2 (left half 8.1 < log(M BH /Mq) < 8.4) and bin 3 
(right half, 8.4 < \og(M BH /M ) < 8.7); (4) log \ E bin 2 (left half -1.0 < log \ E < -0.7) and bin 3 (right half, -0.7 < log \ E < -0.4). 
The symbol on the upper right corner of the each plot shows which sub-bin is considered. The lines are color coded as in Figure [4] 



shapes in different Mbh or logA_e bins show no significant 
difference. 

The lack of SED shape dependency on z, L^i, Mbh 
and X E we observed in the XCRQ407 and SSRQ200 sug- 
gests that neither the emission mechanism nor the accretion 
disk and 'torus' structure alters systematically or dramati- 
cally with these parameters. In other words, a single intrin- 
sic quasar SED in the optical/UV to near-infrared range is 
a meaningful concept. 



4 SED DISPERSION DEPENDENCY ON 
PHYSICAL PARAMETERS 

Even if the mean SEDs show little or no dependency on 
physical parameters, the SED dispersion may change with 
them. We checked the SED dispersion in different bins for 
both the XCRQ407 and SSRQ200. In order to exclude the 
dispersion caused by different brightness of quasars and con- 
centrate on the SED shape dispersion, we only consider the 
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Figure 6. The dispersion of the SEDs for XCRQ407 in bins of z and Lboi after normalization at 1/im compared to E94 radio-quiet SED 
dispersion (red solid curve). 



dispersion of the SEDs after normalizing all the SEDs at 
1/xm. The SED dispersion in different bins at certain spe- 
cific wavelengths are listed as examples in Table [6] 

The resulting dispersion of normalized SEDs in z and 
L bo i bins for XCRQ407 compared to E94 radio-quiet SED 
dispersion are shown in Figure H The XMM-COSMOS 
SEDs generally have a larger dispersion in the optical to UV 
range compared to E94 radio-quiet sample. This is probably 
because E94 is biased toward blue quasars, unlike the XMM 
selected XCRQ407, thus it does not include quasars with 
various optical shapes/colors. The 1-10 /jm dispersion of 
XCRQ407 is closely similar to E94, except at ~ 2fj,m, where 
it is somewhat lower. The XMM-COSMOS SED dispersion 
is almost exactly the same in the 3 /im to 3000 A range for 
different z or Lboi bins. The SED dispersions in this range 
are all below a factor of 2. 

For the host-corrected SSRQ200 sample, the dispersion 
of the SED for each bin is shown in Figure [7] We can see 
that the dispersion in different Mbh bins has most similarity 
between bins compared to other parameters (z, Lboi and 
log As). The SED dispersion in the 3 fim to 3000 A range 
is not as tight as the un-corrected sample, implying extra 
dispersion induced by the host-correction process depend 
on z, Lboi or Xe- The lowest bin for each parameter always 
has the largest dispersion. The 3— 8[im dispersion is larger 
than E94 in the bottom one or two bins for all parameters. 

In all these SED dispersion plots, the UV dispersion 
is generally larger than the near-IR dispersion. This could 
mean that the reprocessing or the hot dust component is 
more uniform than the accretion disk emission. Alterna- 
tively, this could be completely caused by the variability 
of the quasar which affects the UV SED most. 



5 DISCUSSION AND CONCLUSION 

We analyzed the dependence of both the mean and disper- 
sion of the SED shapes in the optical-UV to IR range on 
the parameters z, Lboi, Mbh and logAs for the 407 XMM- 
COSMOS radio-quiet type 1 AGN sample. We also calcu- 
lated the bolometric correction at UV-optical to near-IR for 
the host-corrected SSRQ200, and in four different redshift 
bins. 

As the XMM-COSMOS quasar sample is an X-ray- 
selected sample, it includes also sources with low nucleus 
to host contrast. The mean SED of the whole sample is 
greatly affected by the host galaxy emission for low red- 
shifts and luminosities. Therefore, we mainly studied the 
mean and dispersion SEDs of a host-corrected sub-sample 
of 200 radio-quiet quasars, SSRQ200, in four quartile bins of 
z, Lboi, Mbh and Xe- The mean SED shapes in the differ- 
ent bins are quite similar to each other. These SEDs are also 
generally similar to the E94-like mean SED shape, but they 
tend to have less pronounced optical to UV bump emission 
than in E94, as E94 is biased towards blue quasars. Even if 
we fix one parameter (e.g. redshift), the mean SEDs show 
no obvious evolution with the others. 

We checked the dispersion of SEDs in different param- 
eter bins. The near-IR SED dispersion is generally smaller 
compared to the UV SED dispersion, which might be due 
to the variability of the quasar that mainly affects UV SED. 
The SEDs before and after host correction have very sim- 
ilar dispersion in the 3 /jm to 3000 A range for different 
bins, which implies an invariant intrinsic dispersion of SED 
shapes in this wavelength range. However, in this wavelength 
range, the host-corrected SED dispersion is not as tight as 
the uncorrected sample, probably due to the extra dispersion 
induced by the scaling relationship depends on the physical 
parameters. The host-corrected SED dispersion in different 
Mbh bins has most similarity between bins compared to 
other parameters (z, L bo i and logAe). 





logy(Hz rest frame) logf(Hz rest frame) 



Figure 7. The dispersion of host-subtracted SEDs for SSRQ200 in bins of z, L^oi, Mbh and \e after normalized at l/im compared to 
E94 radio-quiet SED dispersion (red solid curve). 



There is no obvious sign of dependency of the mean SED 
with z, Ltoi, Mbh or logA_g- This result implies that the 
global quasar structures (the accretion disk and the 'torus') 
are invariant with respect to these parameters. Despite the 
strong expectation that feeding of the AGN and host galaxy 
is likely to change (e.g. Merloni et al 2010), as well as the 
X-ray corona is likely to change (aox varies with \e, e.g. 
Lusso et al. 2010), the accretion disk and the 'torus' seem 
to not depend on this mechanism. Thus an intrinsic SED in 
UV-optical to near-IR exists. 

The differences among the SEDs show there might be 
diversity only at the second order. Host contamination, red- 
dening and quasar variability should lead to a dispersion in 
the UV and near-IR SEDs. The lack of SED dispersion evo- 
lution indicates that these effects themselves do not show 
much systematic dependence on physical parameters. 



All the above conclusions depend on the reliability of 
the host correction process. This assumes the scaling rela- 
tionship between black hole mass and bulge mass (e.g. Mar- 
coni & Hunt 2003) with an evolutionary term (Bennert et 
al. 2010, 2011). Uncertainties in this correlation might limit 
our ability to detect SED shape dependency on physical pa- 
rameters. 

A less model dependent analysis will be described in 
a following paper (Hao et al. 2012b), where we introduce 
a new mixing diagram to readily distinguish different SED 
shape, i.e. galaxy-dominated, quasar-dominated and redden- 
ing dominated SEDs. 

So far, we have only studied the SED shape one decade 
on either side of lfim. To study the SED shape evolution 
in UV or FIR is less easy. The UV is greatly affected by 
both reddening and variability that is different for different 
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Table 5. Number of SSRQ200 quasars in Sub-bins for Partial 
Evolution 



bin 
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sub-bin 






sub-bin 


z 
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\ E 




1 
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z(2) 
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13 


14 


16 
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7 


8 
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1 
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10 
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z(S) 
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16 


16 
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lu o ^ooL V / 
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13 




11 


12 




4 


5 




11 


5 




1 


8 




10 
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13 
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12 
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13 
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19 


17 
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16 


11 




15 
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6 
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10 


10 




15 
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1 
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14 
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logA B (2) 
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12 






4 


9 


4 


13 






1 


10 


9 


15 




logA B (3) 
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11 






4 


10 


12 


9 





sources. The FIR is greatly affected by host galaxy star for- 
mation, which also varies from case to case. We will discuss 
these two regions of the quasar SED in following papers. 
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